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From data recently reported from the IceCube telescope, we derive new bounds on the monochro¬ 
matic neutrino signal produced from dark matter particle decays. In the few TeV to tens of TeV 
energy range, these bounds turn out to be better than previous limits by more than an order of 
magnitude. As a result, intensity constraints on neutrino lines at energies above a few TeV are now 
comparable to those on gamma-ray lines. From the same data sample, we also perform a detailed 
search for a neutrino line, showing that there is no significant hint for such a signal. 


I. INTRODUCTION 

The observation of a high-energy spectral line in the 
cosmic gamma-ray spectrum is one of the cleanest indi¬ 
rect detection signatures one could think of for establish¬ 
ing the existence of a dark matter (DM) particle. Such 
a line could be produced from DM particle annihilations 
or decays in our galaxy m- Forthcoming satellite [SHE] 
and Cherenkov telescopes mm, as well as the current 
Fermi Large Area Telescope (LAT) [TT], H.E.S.S. [H], 
MAGIC [13] and VERITAS [fl] instruments, will allow 
to probe this possibility with further sensitivity and to¬ 
wards higher energies. A similar ‘smoking gun’ evidence 
of DM particles would be the observation of monochro¬ 
matic cosmic neutrinos. For various reasons, this possi¬ 
bility has in general been less considered than the one 
of a gamma-ray line. We will here investigate the latest 
IceCube data available to search and constrain neutrino 
line (and line-like) signals from decaying DM. 

On the theoretical side, one could in fact expect many 
decaying DM particle setups to give production rates of 
monochromatic neutrinos that are similar or larger than 
those into photon lines (see, e.g., [T3H25] and for counter 
cases, e.g., pTl 

On the experimental side, however, the traditionally 
reported energy resolution for cosmic neutrinos has been 
far from reaching the level of precision that holds for 
gamma-ray experiments. Still, this argument is not much 
relevant anymore as the neutrino energy resolution now 
reported by the IceCube collaboration is of the order of 
15% for cascade events [30]. It is thus close to what is 
achieved by current gamma-ray experiments in the GeV 
to TeV energy range. Simultaneously, IceCube’s effective 
area has increased, which leads to further improved sen¬ 
sitivity to neutrino fluxes, even after strong data cuts are 
applied to reduce atmospheric backgrounds. Still at the 
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experimental level, one should add that bounds usually 
quoted on the annihilation cross section and decay width 
are much less stringent for neutrino- than for gamma¬ 
lines. For example, in the case of a decay, the most 
stringent IceCube analysis reported so far m quotes a 
lifetime sensitivity of order 10^® sec at TeV energy, or 
few 10^® sec above 10 TeV, to be compared with the 
H.E.S.S. photon sensitivity of few 10^® seconds within 
the 1 — 50 TeV energy range [HI [33] . Given the fact that 
no positive signal has been reported for a gamma-line at 
these energies, such a sensitivity comparison would imply 
that only models which predict much more monochro¬ 
matic neutrinos than photons could be probed by neu¬ 
trino detectors. The main purpose of this work is to 
show that this last statement does not hold anymore for 
decaying DM. 

In our dedicated line search, we will include energy dis¬ 
persions and optimize the statistical method for analyz¬ 
ing the energy spectrum. From recently released IceCube 
data [ 311133 ] , we then determine bounds on the flux of 
monochromatic neutrinos that could have been emitted 
by DM decay in the I TeV to 100 PeV energy range. 
Above 100 TeV DM masses, our bounds are similar to 
the ones recently derived in Refs. [ 131133 ] using IceCube’s 
high-energy dataset from [33]. Below these energies, and 
above few TeV, the bounds we obtain improve previous 
ones by more than one order of magnitude. As a re¬ 
sult, there exists now a region of DM particle masses, in 
the range from several TeV to 50 TeV, where stringent 
constraints on both monochromatic neutrinos and gam¬ 
mas do exist and are comparable. Interestingly, unlike 
at higher energies, this region is relevant for thermal DM 
candidates. 


II. BOUNDS ON MONOCHROMATIC 
NEUTRINOS 

The IceCube collaboration recently released a data set 
of high-energy neutrinos collected during 2 years, from 
2010 to 2012 [32]. At the highest energies, the neutrino 
events partly overlap with a longer exposure - 3 years 
- data set published few months earlier [ 33 ] ■ The new 
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FIG. 1: Deposited energy spectrum from the full sky (black 
curve) as measured by IceCube. The colored regions show 
the expected atmospheric muon (bottom green) and neutrino 
(middle red) fluxes as well as the best-fit astrophysical neu¬ 
trino power-law background model contribution (top blue). 
Data taken from Ref. [32) . 


data set has nevertheless the advantage that it extends 
down to lower energies, which can be of great interest for 
DM searches. From public IceCube data, it is therefore 
now possible to derive the strongest current bounds on 
neutrino signals from DM decays in the TeV to PeV DM 
particle mass range. 

In this study, we will take advantage of the latest 
data of Ref. [52] and carefully derive upper bounds on a 
monochromatic neutrino line signal from decaying DM. 
The bounds we obtain will be compared to the current 
corresponding limits published by the IceCube collabo¬ 
ration [5T] below 100 TeV and to those in Refs. [551155] . 
which instead used the 3-year high-energy data [55] . 
above 30 TeV energies. Other previous bounds can be 
found in Refs. [missii]. 

The energy range covered by the data release in 
Ref. [55] goes from 100 GeV to 10® GeV. Twenty en¬ 
ergy bins have been considered within this range with 
in total iVtot = 383 detected events. The number of 
events in each bin i is denoted by and is shown 

in Fig. [^ We will derive limits on DM signal by com¬ 
bining this neutrino data with the information on the in¬ 
strument response given in the supplementary material 
[55] to Ref. [55]. The search for a possible detection of a 
monochromatic neutrino line and an improved approach 
to derive limits will be done in Section Hi D1 


A. Dark matter signal prediction 

If the DM particle decays into neutrinos, with associ¬ 
ated lifetime it gives rise to a neutrino signal both 


from the Milky Way DM halo and from extragalactic 
DM. The former contributes with an intensity per solid 
angle D and neutrino energy Ey that is 
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The integration is performed along the line of sight, 
where s represents the distance to the observer. This 
intensity depends on the sky direction, which is repre¬ 
sented by the dependence on Galactic coordinates 6 and 
1. However, in the case of a spherically symmetric DM 
halo, the signal intensity only depends on the relative an¬ 
gle 0 between a given sky direction and the halo centre 
direction. We also take the DM halo to coincide with the 
Galactic centre (GG) at (6, 1) = (0,0). The distance from 

the GG is given by r(s, tp) = (rg -I- — 2rs cos ip) and 

we take the Sun to GC distance to be r© = 8.33 kpc 
[451143] . For the DM density we use a NFW profile [44] , 
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with pq = 0.39 GeV/cm® and a = 24 kpc [451147] . 
dN/dEy is the neutrino-energy spectrum at the source, 
which in the case of one monochromatic neutrino per DM 
particle decay would be dN/dE = 5{E — Flune).^ 

As for the extragalactic DM signal, it has an isotropic 
intensity given by 


dEydVl 


47rmDMT'DM 
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where H{z) = -|- Hm(l + is the Hubble ex¬ 

pansion rate as a function of redshift z, c is the speed 
of light in vacuum and pc = 4.769 x 10“® GeV cm“® 
denotes the critical density of the Universe. We fur¬ 
thermore assume a ACDM cosmology with parameters 
Ha = 0.6844, = 0.3156, Hdm = 0.2647 and 

h = i7o/100kms“^ Mpc“^ = 0.6727 — as determined 
from the Planck satellite’s observations and found in Ta¬ 
ble 3 of Ref. [ST]. An estimate of cosmological neu¬ 
trino opacity reveals that emission below redshift z < 
10^(Ai,/TeV)“^/^ undergoes negligible attenuation [55]. 
We therefore neglect absorption effects in Eq. ([^. 


^ Doppler shift and electroweak correction effects broaden 
monochromatic lines. Typically, this broadening is significantly 
less than telescopes’ energy resolutions. Though, for neutrino 
energies orders of magnitude larger than the WjZ boson mass, 
electroweak corrections can be sizable. At leading order, for ex¬ 
ample, Z-strahlung broadens a line so that only 99(89) % of its 
original intensity remains within a 15 % energy band around its 
peak for = 10^(10^) GeV [using in Eq. (11) of |48 | 

for the neutrino line shape, with mv = ^nz and 2 = 2Ejylmjp)M 
in their equation]. For large DM masses, though, higher-order 
resummation techniques might be required [491I50| . We will not 
include these electroweak corrections. 
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Eq. (§ shows that even if DM particles only decay into 
monochromatic neutrinos at the source, the observed cos¬ 
mological signal d4>eg/dE^dU, is not monochromatic due 
to the redshifting effects. Instead of a monochromatic 
line at ifune, it gives a broader line-like spectrum peak¬ 
ing around ifiine/2. For ip > 90° the cosmological DM 
signal intensity is larger than that from the Galactic DM 
halo if integrated over an energy range from ifiine /2 to 
inline- However, at lower latitudes the Galactic contribu¬ 
tion dominates strongly and, e.g., at ip = 10° it is 5 times 
larger than the extragalactic one. As a result, the Galac¬ 
tic signal always dominates when considering the full sky 
flux. For IceCube observations, the Galactic domination 
is further enhanced because IceCube leans towards the 
GC and high-energy (typically E^, > 40 TeV) up-going 
neutrinos undergo non-negligible absorption effects when 
going through Earth. 


B. IceCube’s instrument response 

The differential number of expected events in the Ice- 
Cube detector can be evaluated by convolving the neu¬ 
trino intensities from Eqs. 0 and 0 with the detector’s 
instrument response: 

__ _ d((()h -|- (peg)a p „ / 

dEAndE'dcosB'dcp' ~ dE^dfl “ ^ ’ 

where E', O' and (p' are events’ reconstructed energy. 
Zenith angle and azimuth angle at IceCube, respectively. 
We here also introduce an index a G {e,/r, r, e,/i, f} to 
keep track of the flavor composition of incoming neutrino 
fluxes at Earth’s surface and if they are neutrinos (y) or 
anti-neutrinos {P). 

In general, Eq. Q gives a spatial-dependent energy 
spectrum (induced by the superposition of the spatial 
varying intensity of the Galactic DM signal with the 
isotropic, but cosmologically redshifted, spectrum). 

The exposure Sa = AeB,aT is the product of the ef¬ 
fective area Aeff^a{Eu,0,(p) and the exposure time T = 
J dt = 641 days.^ The Zenith direction z at Ice- 
Cube is, in Galactic coordinates, at bjc = —27.4° and 
lie = 303°, which allows to express the angle rela¬ 
tive to the Zenith by 9{b, 1) = arccos(i • s) for a given 
line of sight unit-vector s{b,l). An incoming neutrino, 
with sky direction {0,(p) and true energy E^ at the sur¬ 
face of Earth, has an expected probability distribution, 
in reconstructed deposited energy E' and reconstructed 
sky direction (9',(p'), given by the dispersion function 

Df^{E' ,9' ,dp'-E,,9(b,lA)AiP,lA))-^ 


A bin integrated version of the differential instrument 
response functions (AeffDeff)ct, pre-integrated over all (p 
and (p', as well as over certain bin ranges in 6, 9' E and E', 
is published in the supplemental material [33] to Ref. [33] . 
This is the instrument response function we use in this 
study.^ It is worth noticing that there is a separate in¬ 
strument response function for each neutrino flavor — 
e, r — as well as if it is v and P, and that they show 
strong dependences on IceCube’s Zenith angle 9ic due 
to energy losses and absorption effects in Earth for high 
energy neutrinos. 

For a given signal in Eq. the expectation value of 
the differential number of observed neutrinos is 


I 11 / __ 
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where we integrate E' over each energy bin AiE' using 
the same binning as in the publication [32]. Because 
Ref. [32] provides no data on (p'{t) and only north/south 
separation in 9', we integrate over all cp' and 9' without 
optimization with respect to the DM signal’s particular 
morphology (i.e. A9'{t) and A(p'{t) include the full sky). 

This leads to our final expression for the expected num¬ 
ber of observed DM neutrino events in each energy bin i: 
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where the incoming neutrino flux’s flavor composition at 
the surface of Earth is set by the probabilities Pa- As 
an example. Fig. [^ shows the energy dispersed signal for 
100 TeV monochromatic neutrinos from DM particle de¬ 
cays (i.e. DAI v-\-X^ where X is any state with negli¬ 
gible mass and with no signal contribution in IceCube). 
In this example, the flavor composition was taken to be 
^e,/i.T = l/3 with equal parts of v and P. 

In general, after propagation, the average flavor com¬ 
position at Earth’s surface is set by the neutrino mixing 
probabilities (using expressions for long baseline oscilla¬ 
tion L y> E^jAmp [S3]): 

P(r'e GG Vf.) = 0.573, P{ve v^) = 0.277 

P{Vf. GG Vt) = 0.150, P{v^ GG = 0.348 

P(v^ GG Vr) = 0.375, P{Vt ^ Er) = 0.475. 


For the above results the following mixing angles have 
been used: sin^ 9i2 = 0.304, sin^ 6*23 = 0.452, sin^ @13 = 
0.0218 [53]. Thus, a pure flavor state (3 produced at 
source gives Pa = PiP ^ a) in Eq. 0 . Depending 
on the flavor of the neutrinos emitted, the bounds on the 
DM lifetime vary by less than a factor two and for a pure 
ly 01 i> flux the limits can vary a bit more — we comment 
more on this in Appendix [C| 


^ We approximate Aeff to be independent on the azimuth angle 
4>{b,l^t) and that the detector has no explicit time dependence. 
^ commonly normalized such that its integration over E', 

cos 9' and (j)' gives 1. 


The effective areas provided in Ref. m are summed over all 
interaction channels. 
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FIG. 2: Deposited-energy (E'dep) spectrum for DM decay into 
monochromatic neutrinos of energy Sjf"® = 10® GeV. The sig¬ 
nal normalization is for a DM particle with mass 2 x 10® GeV 
and lifetime tdm = 10^®s, with democratic neutrino flavor 
mixture and equal parts of v and v in the flux reaching Earth 
surface. The diagram shows the total number of events (black 
curve) after integration over the full sky and including prop¬ 
agation effects through Earth. The dominating Galactic DM 
contribution (blue) and the extragalactic DM contribution 
(red) are shown separately. The dashed curves enclose the ef¬ 
fective area uncertainty range when the presented uncertain¬ 
ties of IceGube’s instrument response function in the different 
sky directions are propagated and added in quadrature. The 
energy binning is the same as in the published instrument 
response function |32| . 

C. Robust dark matter constraints 

From the observed number of neutrino events in each 
energy bin Ni (given in Fig.[^, robust upper bounds can 
be put on any model’s neutrino flux prediction. 

At the <7 = 95 % confidence level (CL), an upper bound 
AYjjnit in each bin is, by Neyman’s construction [5S] from 
a Poisson distribution, given by® 

fc =0 

There is currently no certain understanding of the as- 
trophysical neutrino background contribution at these 
energies. Therefore, we will only include the expected 
muon induced and atmospheric neutrino backgrounds. 
These background were derived in Ref. m) and the 


^ We note that our classical frequentist approach differs from the 
Bayesian approach in Ref. m (and potentially also from m), 
but numerically the results are practically identical except at 1-2 
TeV DM masses where their Bayesian approach would give limits 
marginally weaker by 30%. 



FIG. 3: 95 % CL limits on the lifetime of DM particle decay 
into monochromatic neutrinos. Solid (black, blue and red) 
curves are for DM decay into pure neutrino flavor states {ve, 
Vfi and v-r, respectively) at production - which after propaga¬ 
tion to Earth are no longer pure flavor states. The companion 
final state X in the DM decay DM Va + X is assumed 
to be light compared to the DM particle mass moM- The 
dashed red curve shows the limit for DM —>■ i/,- -b A if no 
atmospheric background subtraction is made. For compari¬ 
son, we show the limits by Rott, Kohri and Park (90 % 
CL, green dotted curve) as well as Esmaili, Kang and Serpico 
[34j (90 % CL, green dashed-dotted curve) who analyzed the 
3-years high-energy data set from Ref. |35| . Equal parts of u 
and u are assumed for the DM signals. 

sum of their central values and respectively, 
are represented by the upper edge of the red region in 
Fig.[T} To derive robust limits, we scale down these 
backgrounds to their 2-sigma lower edges, such that 
Af^kg = 0.538Af^ + 0.938A®. The 1-sigma uncertainties 
of the penetrating muon N'^ and atmospheric A® back¬ 
grounds are read from Table 1 in Ref. [32] • The lower 
bound Tiiniit on the DM lifetime is then set so that the 
DM plus background, expectation values 

do not overshoot the experimental limits AYimit io any 
energy bin. For each DM mass we thus set 

Rimit = iiiin{TDM € : Afpjyj -|- A^kg < Al]®;jj^k.}. (8) 

Note that A^kg < Af®ii„it always fulfilled in this data 
set.® 


® A violation of this inequality would be interpreted as an arti- 
fact of setting the background to a fixed assumed known value — 
whereas it should be associated with a certain uncertainty - or 
as a large statistical downward fluctuation in the data. Statisti¬ 
cally, one could also worry that the limits are derived from the 
most constraining bin. However, in practice it is at most two E' 
bins that are relevant and a trial factor of two would at most 
make a 95 % CL drop to a 1 — (95%)^ ~ 90 % CL limit. 
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In Fig. we give 95 % CL bounds on DM decay into 
a monochromatic neutrino line. We show limits for the 
cases of pure v,., or states and with equal parts of 
anti-neutrinos produced by the DM decay, DM — Va+X. 
For illustration, we also show the limits if zero back¬ 
ground is assumed (thin dashed curve), although this is 
unrealistically conservative. In this hgure, the bounds 
obtained are compared with the ones from Refs. [231 
[31j, in the range from 60 TeV to 10 PeV DM masses. For 
the DM signal, we assume equal induced fluxes of neutri¬ 
nos and anti-neutrinos. The peak at muM = 12-6 PeV 
seen in the figure arises from the Glashow resonance of 
6.3 PeV electron-antineutrinos annihilating on electrons 
in the ice [56ll57] . 


D. Search for a i^-line signal and improved bounds 

1. Method 


Besides setting the best possible bounds on the DM 
lifetime, it is also of interest to search for any hint of a 
line signal in the data. We propose a novel methodology 
for neutrino line(-like) searches: a profile log-likelihood 
study of the energy spectrum including proper energy 
dispersions. This differs from previous studies that used 
Kolmogorov-Smirnov tests or approaches with no de¬ 
tailed spectral shape analysis na m uniM iMsmsH] . 
The method is similar to, e.g., the Fermi-LAT gamma- 
ray line search [TT] (also briefly summarized in [33] )■ This 
approach will enable us to search for line (and line-like) 
signals and further improve the DM bounds of the pre¬ 
vious section. 

To claim a detection of a monochromatic line requires 
enough statistics and energy resolution to distinguish 
such a feature in the spectrum from an assumed back¬ 
ground model. A simple model to fit the neutrino flux 
that exceeds the expected muon and atmospheric back¬ 
grounds is an isotropic astrophysical signal with a single 
power-law energy spectrum, with equal parts of each fla¬ 
vor as well as of neutrinos and anti-neutrinos, 


d<^astro 


= 3 X 10 "^® X ^0 



( 9 ) 


where Eq = 10® GeV. Such a power-law spectrum can 
be expected from conventional astrophysics contribu¬ 
tions such as, e.g., active galactic nuclei and start burst 
galaxies with 2.0 < 7 < 2.6 [301 ISI|- The best-ht val¬ 
ues from [32] are 7 = 2.46 ± 0.12 and (po = 2.06toJ6 
GeV“^cm“^sr“^s“^. To quantify the goodness-of-fit, we 
performed a Pearson test [32], 
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bins i 


i^ohs Pi^tot)'^ 

tot 
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Here, the set of Pi = Nl^/Ntot gives the deposit-energy 
probability distribution of the model. The NIq = -|- 


Nl + are represented by the upper edge of the 

blue region in Fig. [^ This will be our NULL model for 
generating Monte Carlo data. The IceCube data gives 
X^ = 12.7 for this NULL model when we include all the 
20 energy bins from Ref. [32]. Pseudo experiments on 
10® Monte Carlo-generated data sets show that this cor¬ 
responds to a P-value of 0.42.^ The conclusion is that 
this models give a good fit and there is no apparent need 
for a contribution from a DM signal in the current Ice- 
Cube spectrum. 

Given that such a simple background model fits the 
data well, we construct a test statistic (TS) by comparing 
the maximum log-likelihood of this type of background 
model to the hypothesis of a monochromatic neutrino¬ 
line signal on top of the background. Hence, to study if 
a line signal improves the ht, we evaluate 


TS = 21n 


■^(u-sig — U-sig,best) 
/l(»^sig = 0) 


( 11 ) 


where the Poisson likelihood function is 


^=n 

bins i 


(v. 


model/ _ 

KJ- 


( 12 ) 


and the flux model is the superposition of the, 

detector convolved, signal and backgrounds contributions 


lV^odel( *T-sig) '^1,2,3) 7 ) ''T'sig-^DM(^DM) T'o) 

+ niN^^-\-n2Nl + , (j)o)■ ( 13 ) 

The three backgrounds’ normalizations rii 2,3 and the as¬ 
trophysical power-law spectral index 7 are set to their 
best-fit (largest likelihood) values for each value of a DM 
neutrino signal amplitude nsig. The normalizations of 
the DM Vpjyj and background signals are relative 

to their (arbitrary) reference points of tq = 10^® s in 
Eq. Q and (j)Q = 2.06 GeV“^cm“^sr“^s“^ in Eq. 
respectively. We note that our best-fit (i.e. maximum 
likelihood) background-alone model has ni = 1.36 ,712 = 
0.85,713 = 1.12 and 7 = 2.62 in Eq. (13). This back¬ 
ground model is plotted in Fig. [^ of Appendix [^ after 
convolution with the IceCube’s instrument response. 

We use the full energy range of data for each TS de¬ 
termination. Because the statistics is fairly limited in 


^ A direct assumption of a distribution, with — 5 = 15 
d.o.f., for the quantity in Eq. ( |10[ | would instead have given a 
P-value of 0.63 for the NULL model. However, because of low 
statistics this assumption is not fully valid, and we generated 
data realizations from Poisson distributions of the number of 
events in each deposited-energy bin to calculate P-values. All 
fits to Monte Carlo data were done with an adjustable single 
power-law astrophysical component 7 ) together with 

free normalization of the two atmospheric background compo¬ 
nents and A' - as in Eq. ( |13| l . In Ref. |32| the goodness-of-fit 
P-value is instead found to be 0.2 for this model; but they con¬ 
sider three observables (not only reconstructed deposited energy, 
but also directional and track property information) in their fits. 
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FIG. 4: TS distribution from 2.02 x 10^ monochromatic neu¬ 
trino line searches (at our 101 used DM masses) in Monte 
Carlo generated pseudo data from our NULL background 
model assumption. The data follows well a |5(TS)-I- |x^(TS) 
distribution (shown by the dashed blue line after multiplica¬ 
tion by 2.02 X 10^). Dotted grey curves show the TS distribu¬ 
tion from moM < 560 TeV (dark grey) and ttidm > 560 TeV 
(light grey) line searches separately. 



FIG. 5: Local fit significance of a monochromatic line from 
DM decays versus DM mass (solid red curve). For comparison 
we also show the TS evaluated for signals at our ‘robust’ limits 
with the method of Section [ll C| (black dashed) and our used 
profile Log-likelihood ratio limit at TS = 2.71 (blue dotted 
curve). Democratic flavor composition and equal parts of u 
and V are assumed at Earth’s surface. 


all the 20 energy bins, we do not expect the systematic 
uncertainties to severely dominate over the large statis¬ 
tical uncertainties even if a single power-law model is 
assumed to be appropriate for the astrophysical back¬ 
ground in this large energy range. By Wilks theorem 
[55] (or rather Chernoff theorem [Bl], as we constrain all 
signal contributors to have non-negative normalizations), 
the TS should asymptotically follow a |d(TS) -I- |x^(TS) 
distribution. However, due to the low statistics, this is 
not guaranteed and we performed 2.02 x 10^ Monte Carlo 
pseudo experiments. For each DM mass range tested, 
we find good agreement between the obtained TS dis¬ 
tribution and the asymptotically expected distribution. 
Figure [^illustrates this by showing the collective TS dis¬ 
tribution from our Monte Carlo data realizations of the 
background NULL hypothesis. 


2. Line search and limit results 

No significant line signal was found in the IceCube 
data. We tested DM masses from 2 TeV to 100 PeV 
in steps of 5 % in logj^Q Enne and checked the various fla¬ 
vors and neutrino/anti-neutrino compositions presented 
in the figures below and in Appendix 0 The maximal 
TS was 2.9, corresponding to a Gaussian equivalence of 
~ 1.7 a significance, and occurred for a line of 
pure r'e at Earth and DM mass of 45 TeV. In Fig. we 
show a representative plot of the significances (solid red 
curve) as a function of DM mass. The plot is for the 
case of democratic flavor composition and equal parts of 


v and P at the DM source. In this setup, the maximal 
significance is a/TC ~ 1.5 a for ttidm = 44.8 TeV when 
tdm = 1.0 X 10^® s. We note that the second largest 
TS peak, with \/TS ~ 1.3 u (for ttidm = 2.52 PeV with 
'''DM = 3.9 X 10^® s), is due to the few neutrino events at 
~ 1 PeV energies for which a DM origin has been specu¬ 
lated |2lH25l|3ll4ll|65H2Q]. 

Our best-fit DM spectrum as well as our best ~ 1 PeV 
line signal are shown for illustration in Fig. j^of Ap¬ 
pendix]^ The 3-year IceCube data sample in (35| and a 
DM signal with an accompanying continuum of lower en¬ 
ergy neutrinos to a monochromatic neutrino line should 
not significantly alter these significances [21], unless the 
background modelling is very different. 


With no significant line detected, we turn to derive 
stringent upper limits on DM decay into a monochro¬ 
matic neutrino line. We set our 95 % CL lower limits on 
t’dM) by requiring that the profile likelihood £(nsig) with 
respect to its maximum values to have TS < 2.71 from 
Eq. ([n]).® The derived limits are shown in Fig. [^ where 
we also indicate by a thin red line the 99.9 % CL limit. It 
turns out that this method improves lifetime constraints 
by up to a factor 5 compared to the approach of Sec¬ 
tion IIC In fact, the shown 99.9 % CL limit follows fairly 
closely the ‘robust’ limits from Fig. [^ All these limits 
are for the case of democratic neutrino flavor composi¬ 
tion Pa = 1/3 and equal part of v and P (see however 


In Appendix[^we verified that this method has the statistically 
correct 95 % coverage. 
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FIG. 6: 95 % CL lifetime limits (solid curve) on the DM par¬ 
ticle decay lifetime into monochromatic neutrino. Expected 
sensitivity reach (dashed curve) and its 68% (yellow) and 
95 % (green) containment bands are also shown. The thin 
red line shows shows the 99.9 CL limits. Democratic flavor 
composition and equal parts of v and h are assumed for the 
incoming neutrino flux at Earth surface. 



FIG. 7: 95% CL lifetime limits on various neutrino line¬ 
like signals from DM decay: monochromatic line, internal 
Bremsstrahlung dN/dE^ oc Ei, and on E^) and box-like spec¬ 
trum {dN/dEi, = constant). Democratic flavor composition 
and equal parts of u and v are assumed after propagation to 
Earth. 


Appendix for variations). 

In Fig. we also compare the derived limits with ex¬ 
pected sensitivity. For data with low statistics, large fluc¬ 
tuations can naturally appear between different data re¬ 
alizations and this comparison becomes particularly rel¬ 
evant. 200 Monte Carlo data sets are generated from the 
NULL hypothesis at each DM mass. The expected 68 % 
(yellow band) and 95 % (green band) limit containment 
bands are shown together with the median expectation 
(dashed black curve). 

For illustration, we calculated the TS at the ‘robust’ 
limit from Section IIC and show them in Fig. by the 
dashed black curve (in this case we use democratic flavor 
composition of the DM signal, but otherwise the limit 
is as in Fig. [^. Except at the lowest masses, they al¬ 
ways have TS values larger than 2.71 (which justifies 
to call them ‘robust’). The reason that they go below 
■\/TS = V2.71 = 1.65 at the lowest DM masses is an 
effect from that our ‘robust’ limit derivations assume a 
fixed atmospheric background contribution while the pro¬ 
file log-likelihood setup allows free normalization of all 
backgrounds. At the lowest line energies, it happened 
that our best-fits prefer lower atmospheric background 
contribution compared to the fixed NULL model from 


3. Constraints on other line-like signals 

The method described in Section al D l| is fully general. 
In particular, many line-like spectra can pragmatically be 


parameterized as a single power-law with exponent n: 


dN 

dE 


— ) 0(mDM-2£;), (14) 

WdM V™DM/ 



where the sharp spectral cut-off at if = todm/ 2 is param¬ 
eterized by the unit-step function 0 (todm ~ 2E). 


This spectral shape works as a good proxy for several 
physical scenarios. It appears in the cases of DM de¬ 
cay into three-body final states from effective operators 
(where n « 2—3 [7T] ) or when DM particles decay into an 
intermediate state that then promptly decays into, e.g., 
neutrinos to form a ‘box-shaped’ spectrum (where n = 0 
[72]). For gamma-ray DM signals, final state radiation 
and internal bremsstrahlung processes also give spectra 
approximately of this form (with n « 0—3 [75H75] ). Elec- 
troweak corrections to a monochromatic neutrino line, 
which we commented on in footnote [l] is however not of 
this form. Instead, the neutrino spectrum gets a low en¬ 
ergy bump with a very steep rise towards its peak energy. 
Nonetheless, for all practical purposes this is equivalent 
to a monochromatic signal. 


The spectra in Eq. (14) are normalized to give one 
neutrino per decaying DM particle. However, e.g., the 
box-shape spectrum is typically thought of as two inter¬ 
mediate states that promptly decay into two particles (in 
this case, two neutrinos), hence four neutrinos per DM 
particle decay with maximum energy of mDivi/d might be 
more natural in this case. We derive limits for the cases 
of n = 0,1,3 for these types of DM signals and show 
those, together with the pure monochromatic line signal 
limit, in Fig. 

















4- Possibilities of improvement 

In this study, we have presented a methodology to 
search and derive limits on a DM-induced neutrino line 
signal. We have always considered the sum of all events 
(e.g. track- plus cascade-like events), the full sky region 
and events binned in energy. Future improvement might 
be possible by using additional information in the data: 

• ‘Morphology sorting’ of the neutrino events, e.g, 
into cascade- and track-like events and correlate 
this to the expected DM and background signals 
[76]. Typically, a muon neutrino produces a track¬ 
like signature, while tau- and electron-neutrino in¬ 
teractions have cascade-like morphologies |77| . Uti¬ 
lizing the method of Section [ll C[ and assuming no 
backgrounds, limits from track and cascade events 
separately did not significantly improve limits. 

• Using directionality information of incoming neu¬ 
trinos to, e.g., find an optimized sky regions where 
signal to background ratio is expected to be the 
largest m [38] . Unfortunately, the publicly pro¬ 
vided events and the effective areas have limited 
informations on this [33] . A simple check with sep¬ 
aration into north and south hemisphere subsets 
did not reveal significantly improved limits. 

• Exploring the data unbinned in reconstructed ener¬ 
gies. This would require a better resolution of the 
instrument response function than the limited 20 
energy bins provided between 10^ to 10® GeV [33]. 

Improvements in sensitivity could therefore be possible, 
but care must be taken to not reduce the already low 
statistics in any analysis. Systematic effects, e.g. in the 
effective areas and the background modelling, might also 
be studied in more detail ISI1I33I35]- A straightforward 
propagation of effective areas uncertainties, taken from 
|33j . by adding them in quadrature revealed less than 
30 % rescaling of limits. 

If we recall that this analysis is based on a data set 
collected only from year 2010 to 2012, the method ap¬ 
proach presented here could even open up the potential 
to discover a DM induced neutrino line signal in already 
collected data with the used, or improved, event selection 
cuts. 


III. CONCLUSIONS 

We have presented in Figs. andnew bounds on the 
lifetime of a DM particle decaying into monochromatic 
neutrinos or similar line-like signals. In Appendix 
we also show complementary limits, including various as¬ 
sumptions on the neutrino flux composition of the emit¬ 
ted signal. Above few TeV DM masses, our derived limits 
are significantly stronger than the latest limits from the 
IceCube collaboration |3T], which we illustrate in Fig. 
For example, around 20 TeV DM masses our limits are 



FIG. 8: Lifetime limits on DM particle decay into monochro¬ 
matic lines V + y. From Fermi-LAT [78] (blue, dashed) and 
H.E.S.S. [121128) (blue, dotted) using gamma-ray data com¬ 
pared to the neutrino line bounds derived in this study (solid 
black) as well as previous IceCube limits |31) (thick green). 

a factor 25 better than the IceCube collaboration’s 90 % 
CL limits (that used a slightly different DM density pro¬ 
file) [31]. For DM masses above 100 TeV and up to 
1 PeV, our bounds are in agreement with those obtained 
in Refs. [531 El] (which were compared to our ‘robust’ 
limits in Fig. [^. The improvements in limits are a con¬ 
sequence of analyzing the latest data published by the 
IceCube collaboration [551 [33] and a refined approach to 
neutrino line searches by studying the energy spectrum 
with a profiled likelihood method. 

For energies above a few TeV, the sensitivities reached 
on neutrino lines are now comparable to those existing 
on monochromatic gamma-ray lines. As illustrated in 
Fig.i the following picture emerges: 

i) For DM masses below few TeV, constraints on decay 
lifetimes to gamma lines remain orders of magnitudes 
stronger than those to neutrino lines; 

ii) Above the maximum energy considered by the 
H.E.S.S. experiment, = 25 TeV, there are to 
our knowledge no numerically precise gamma-line 
constraints (see however [ilTDllTM]) but strong 
neutrino-line constraints exist now up to energies 
several orders of magnitude higher; 

iii) In the multi-TeV to 50 TeV mass range, the lifetime 
constraints for these two monochromatic decay chan¬ 
nels only differ by a factor of one up to an order of 
magnitude. 

With foreseen improvements in both neutrino [53H57] and 
gamma-ray [881191] data, this opens up increased chances 
to see a ‘double-barreled smoking gun’ signal in the form 
of monochromatic neutrino line plus gamma-ray line from 
DM particles m- 
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Appendix A: Best-fit spectra 


To illustrate how some of the best-ht deposit energy 
spectra compare to the IceCube data, we show in Fig|9|i) 
our best-ht background model [presented below Eq. (13)] 
ii) our best-ht DM signal model and hi) our best-ht DM 
model with a monochromatic neutrino line signal around 
the 2 observed events at ^ 1 PeV energy. These are 
spectra in the case of democratic havor and equal parts 
of V and i>. 

The ~ 1 PeV DM signal might seem too low to give the 
best-ht, but we checked that an increased signal would 


very slightly worsen the likelihood of the ht. We also 
show the IceCube collaborations [32] best-ht model from 
Fig0 Similarly to their best-ht background model (con¬ 
sidered in Section IID 1) [32] , our best ht background 
model gives a = 12.4 from Eq. (101. This corresponds 
to a P-value of 0.40 from a Monte Carlo validation (as 
opposed to a P-value of 0.65 if a distribution were 
assumed - c.f. footnote [^. 


Appendix B: Statistical coverage 


The statistical coverage of our limits was investigated 
by simulating a DM line signal on top of our NULL 
hypothesis (given in the beginning of Section IID 11. 
We simulated 10^ Monte Carlo realizations, and looked 
at 101 DM masses between 2 TeV and 200 PeV. The 
monochromatic neutrino lines were randomly given a sig¬ 
nal strength corresponding to a rpM between one order of 

magnitude larger or smaller than our derived DM limits 
_ic 

^ limit ■ 

For lifetimes shorter than our derived limits from the 
IceCube data (^DM/Tumit < 1)^ ^^e found coverage is 
93%, which is in good agreement with our stated 95% 
CL. For TDM/T/imit > 1 the coverage is 99%, which cor¬ 
responds to a safe overcoverage that is expected for low 
signal strengths. In Fig.[^we show the hrst 1010 points, 



FIG. 9: Observed deposit energy spectra compared to: our 
best-ht model with a DM particle mass of 44.8 TeV (blue, and 
blue-dotted for the DM contribution); our secondly best-ht 
model with a DM particle mass of 2.52 PeV (red, and red- 
dotted for the DM contribution); our best-ht background-only 
model (orange); and the IceCube collaboration’s |32| best-ht 
background model (grey). The best ht DM decay signal con¬ 
tributions are shown for the 44.8 TeV and 2.52 PeV DM par¬ 
ticles by the dashed red and dashed blue curves, respectively. 
DM models are for democratic havor composition and with 
equal parts of v and v hux at Earth surface. The signihcance 
for both the DM signals are less than 1.5 ct. 



/ IC 

rpM /T limit 


FIG. 10: Comparison of the 95% CL lower limit 
to the true DM lifetime value tdm for a set of 100 Monte 
Carlo realizations (for each 101 DM masses tested between 
2 TeV and 200 PeV). Each realization is from a model with a 
DM particle lifetime tdm randomly drawn between 0.1 and 10 
times the lifetime limit derived from the IceCube data 
The color of each point gives the found signihcance UTS (i.e. 
the number of standard deviations a) for the injected DM 
model. Points above the solid horizontal line represent real¬ 
izations where the lower lifetime limit covers the true injected 
value. Democratic havor and equal parts of i' and P compo¬ 
sition was assumed at Earth surface for the DM signals. 
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FIG. 11: 95 % CL limits on the lifetime of DM particle decay 
into monochromatic nentrinos for pure e (black), /r (blue) or 
r (red) neutrinos injected at Source and with equal parts of 
1 / and v. 


FIG. 12: Same as Fig. with the pure Ve, i'll and Vr 
at Source, but also separated into neutrinos (solid) or anti¬ 
neutrinos (dashed) 


where we color-code each simulated DM signal with its 
corresponding TS value. For a clearer color scale, we 
assigned points with \/TS > 6 the value 6. 


Appendix C-. v ws u and flavor composition 
dependences on DM lifetime limits 


In this appendix we collect some complementary limits 
that might be of interest. All limits in this Appendix are 
derived by the method described in Section [IlD 1| 

Figure 11 shows the 95% CL lower limits on DM par¬ 
ticle lifetime into a monochromatic neutrino in the cases 
of pure electron, muon or tau neutrinos injected at the 
DM source and with equal parts of v and F. The neu¬ 
trino flavor oscillations from Section Pl Bl are used for the 
propagation from source to Earth’s surface. 


In Fig. [T^ we use the same setup, but assume only 
neutrinos or anti-neutrinos are injected by the DM decay. 


For Fig. 13 we assume that the incoming fluxes at 
Earth are 100% pure Vfi or Ur. As in the previous 
figure, we show limits separately for neutrinos and anti¬ 
neutrino injected at DM decays. These (hypothetical) 
pure flavor fluxes incoming to Earth could potentially be 
partly mimicked by non standard flavor oscillation effects 
from: neutrino decay [HU US], sterile neutrinos [M], a 
pseudo-Dirac nature of neutrinos |95j , Lorentz or CPT vi¬ 
olations [96] or gravity-induced decoherence m- In any 
case, these limits clearly illustrate the main differences 
in sensitivity reach with this data set for extreme cases 
(i.e. for pure flavor and particle/anti-particle fluxes). Po¬ 
tentially, these limits could be further improved if flavor 



FIG. 13: Same as Fig. [^but in the case of pure r'e, and 
Vr fluxes at Earth’s surface. A pure flavor flux is typically 
not realistic due to neutrino oscillations, but the plot clearly 
illustrates the differences in effective area depending on the 
incoming neutrino flavor. 


and anti-neutrino identification properties in the data 
was used (as is mentioned in Section IID 41. 


Appendix D: Tabulated limits 

Tables of all limits and TS values in this article 
can be found in the online supplementary material file 
Appendix_D_Table. dat [included among the source files 
submitted to the arXiv at http://arxiv.org . 
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